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f^ ■ Abstract 

We consider some bilinear recurrences that have apphcations in 
number theory. The expUcit solution of a general three-term bilinear 
recurrence relation of fourth order is given in terms of the Weierstrass 
sigma function for an associated elliptic curve. The recurrences can 

QQ ' generate integer sequences, including the Somos 4 sequence and ellip- 

CO , tic divisibility sequences. An interpretation via the theory of integrable 

||0 ' systems suggests the relation between certain higher order recurrences 

and hyperelliptic curves of higher genus. Analogous sequences associ- 

C^^ ' ated with a q-discretc Painleve I equation are briefly considered. 
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1 Introduction 

Question: What's the next number in this sequence: 

1,1,2,3,...? 

If you thought the answer is 5, then you guessed wrong! 
Actually the required answer is 7, and it continues 

1,1,2,3,7,23,59,314,..., (1.1) 
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but it's a terrible question, as any four numbers can serve 
as initial data for a recurrence of fourth order or higher. 
Unfortunately this is the sort of question that is routinely 
given to schoolchildren - 1 noticed it in a recent GCSE exam 
set to 15-year-olds in Britain. 

So if it's not Fibonacci, which comes from a second order 
linear recurrence relation, then how is the sequence ( 11.11 ) 
generated? It comes from the bilinear fourth order recur- 
rence 

Tn+2 = , To = Ti = r2 = Ts = 1, (1.2) 

Tn-2 

and goes by the name of the Somos 4 sequence in Sloane's 
catalogue P3]. It turns out that with the special choice 
of initial conditions as in ( 11.21 ). all numerators cancel with 
their denominators in such a way that all the terms in the 
sequence are integers. 

The remarkable fact is that the sequence ( 11.11 ) is related 
to the elliptic curve 

y^ = 4x^-4x+l. (1.3) 

Morever, there is a whole family of sequences defined by 
fourth order bilinear recurrences, of the form 

rn+2rn-2 = « Tn+iTn-l + /? (t^)^ (1-4) 

with constant coefficients a, p. These include a type of 
integer sequences, known as elliptic divisibility sequences, 
P!B1 [TBj , and their generahzations |T2l [15] • The connection 
with elliptic curves is a general property of the family of 
fourth order bilinear recurrences. The following theorem is 
proved in [7j. 

Theorem. The general solution of the quadratic recurrence 
relation ( [i.^[ ) for a ^ takes the form 

r„ = AB^^-^±^, (1.5) 



cr(/«)" 



where k, and zq are non-zero complex numbers, the con- 
stants A and B are given by 

A = ^,, B = ^JfM2>}li, (1.6) 

a{zo) a{zo + k) tq 

and o denotes the Weierstrass sigma function associated to 
the elliptic curve y^ = ix^ — g2X — g^. 
The complex constants zq, k should be considered as points 
on the Jacobian of the curve, hence only defined modulo the 
period lattice, and together with the invariants g2: gs they 
are determined uniquely in terms of the constants a, P and 
the initial data tq ^ 0, ri, r2, ra specified for the fourth order 
recurrence (11.41 ) . The precise construction of the curve from 
a, (3 and the initial data is presented in [7]. Note that in 
the case a = excluded by the Theorem, the recurrence 
( 11.41 ) decouples into two trivial recurrences for even and 
odd terms, so that T2k = tq{t2/ tq)^ (3^^^^^^ ^"^ ^ and T2k+i = 

The terms in the Somos 4 sequence ( 11.11 ) are given by 
formula ( 11.51 ) with 

/oo 
(4^3 - 4t + 1)"^ dt = 1.859185431 

/oo 
{4:t^-4:t+l)--^ dt = 0.204680500+1.225694691Z, 

where the integrals are taken along suitable paths in the 
t plane avoiding the three real roots of the cubic, and 
have been evaluated to 9 decimal places using version 8 
of the MAPLE computer algebra package. The real and 
imaginary half-periods are cji = 1.496729323 and cu^ = 
1.225694691^ respectively. The other constants appearing 
in ( 11.51 ) are just 

a{K) = 1.555836426, A = 0.112724016 - 0.824911686i, 



B = 0.215971963 + 0.616028193^. 

However, in this case the sequence of arguments of the 
sigma function can be written more succinctly as 

ZQ^nK= (2n - 3)io, zq = 0.929592715 + cja, 

so that the iterates of the recurrence correspond to the 
sequence of points (2n — 3)P on the curve (11.31 ). where 
P = (0, 1); a referee was helpful in pointing this out. 

There is currently much number theoretical interest in 
elliptic divisibility sequences (integer sequences satisfying 
Tn\Tm whcncvcr n|m), and in Somos sequences, since ap- 
parently they are a source of large prime numbers [H, [5] . 

2 Symplectic structure and integrability 

By the change of variables 

In- (^__j2 . ('■') 

the bilinear recurrence is related to a second order differ- 
ence equation, and so Tn is a tau-function. 
Proposition. For a 7^ the second order nonlinear map 

Jn—ljn \ JnJ 

has the general solution 

fn = p{k) - p{zo + tik) (2.3) 

with complex parameters zq and k ^ 0, where the p is the 
Weierstrass elliptic function associated to the curve 

y^ = 4x^ - g2X - 5^3, (2.4) 



and the points (p(/t), p'{k)) and {p{zo), p'(^o)) on the curve, 
together with the values of the invariants g2, gs, are deter- 
mined explicitly by the parameters a, (3 and the non-zero 
initial data /o, /i of the map ^KE)- 

The second order discrete equation (12.21 ) can be thought 
of as a mapping of the plane, {fn-h fn) -^ (/n, /n+i), with 
the first integral 
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satisfying X{fn, fn+i) = Kfn-i,fn)- This two-dimensional 
mapping is not symplectic, because {fn+i)~^dfn A dfn+i = 

{fn-iy^ dfn-l A dfn. 

To get a symplectic map, we make the transformation to 
symplectic coordinates {q^Pn) on the curve (with two-form 
dpn A dqn)^ by setting /i = y/a and 

9n = A - fn, Pn = -fnifn+1 " fn-l)- (2.5) 

/i 

With these variables the mapping can be rewritten as 

Qn+l = ^„ - 9n - A, .^ g^ 

Pn+l = -2^1 + {Mn + 2A)$^ - Pn, ^ ' ^ 

where 

^n = ^/ _ ^y KQn^Pn) = V PI + ^(A^ - Q^) - 5^2 (A - Qn) ■ 

Note that A and g2 should be interpreted as constant pa- 
rameters for the map (12.61 ). while /i is interpreted as a func- 
tion on phase space. 

The symplectic map (12.61 ) is integrable, with first integral 



gsiqn^Pn) = Ml - g2qn-pl, (2.7) 



and has a 2 x 2 discrete Lax pair [7J . With a different choice 
of symplectic structure on the eUiptic curve, it is equivalent 
to the one-particle discrete Garnier system in jB]. 

It is natural to wonder how all of the above might gen- 
eralize to the case of the bilinear recurrences 

9 

Tn+g+lTn-g-l = / ^ Oij Tn+jTn-ji (2-8) 

which for 5^ > 2 is the natural higher order analogue of 
( 11.41) . with g+1 parameters ao, Q^i, • • • , %• The arithmetic 
of sequences defined by such recurrences is discussed in [B] . 
We are able to state the following 

Conjecture. The general solution of the bilinear re- 

currence 112. S[) for g > 2 is given in terms of the Riemann 
theta-function associated to an algebraic curve F of genus g. 
Using the change of variables ^2.1\) . the bilinear recurrence 
yields a nonlinear discrete equation^ of order 2g, which is 
equivalent to an integrable symplectic map that linearizes 
on Jac(r), the Jacobian of the algebraic curve. 
The complete proof of this conjecture is work in progress 
with Harry Braden, but here we briefly outline its main in- 
gredients and explain why the result is plausible. The result 
can be seen as a direct consequence of certain vector addi- 
tion theorems for theta-functions obtained by Buchstaber 
and Krichever [1]. However, it turns out that the recur- 
rence ( 12.81 ) is naturally related to a hyperelliptic curve F of 

the form 

2.9 

^2 = 4^25+1 ^^^^.^2g-,^ (2.9) 

(the KdV spectral curves [S]) rather than to an arbitrary 
curve of genus g. Furthermore, it is most convenient to ex- 



The recurrence is ft^ YlltlUn+kfu-kY^^-" = ao + 



Ef=l«./^nLl(/n + fc/n-fc)^-'=. 
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press the general solution using the Kleinian sigma function 
of the curve (this is equivalent to an expression in theta- 
functions - see [2] ) • In that case the formula for r„ is almost 
identical to ( 11.51) except that zq G Jac(r) is a vector in the 
gr-dimensional Jacobian corresponding to the reduced divi- 
sor Dq = (Pi — oo) + (P2 — 00) + . . . + (Pg — 00) of g points 
on the curve, and k G Jac(r) corresponds to the reduced 
divisor D = P — 00 of. another point P (with a derivative 
of a taken on the theta divisor). 

The iterates of the recurrence correspond to the sequence 
of divisors Dq + nD, which gives a linear flow zq + uk on 
Jac(r). The division polynomials for hyperelliptic curves, 
corresponding to the special case of multiples of a single 
point, nD = n{P — cx)), have been considered by Cantor 
[3], who gave examples of these recurrences. Matsutani has 
also done some detailed calculations of addition formulae 
in genus two [ID] . 



3 An integrable g-difFerence analogue 

The second order nonlinear difference equation (12.21) has a 
non- autonomous, g-difference version given by 

Jn+lJn-l — p • i<J--L) 

J n 

This map satisfies the singularity confinement test and has 
a 4 X 4 Lax pair jTTj, as well as a continuum limit to the 



first Painleve equation (PI). Therefore (13.11) may be referred 
to as qdPI. The bilinearization of qdPI is achieved by the 
same formula (12.11 ). to obtain the the g-difference analogue 
of (O), that is 



n _ _ , /Q /_ n2 



rn+2rn-2 = a^" Tn+iTn-i + /? (r„) . (3.2) 



It is a remarkable fact that the bihnear equation ( 13.21) 



shares some of the divisibility properties of the autonomous 
recurrences ( 11.41 ). in the sense that given constants r_2, r_i, 
To, Ti as initial data, all of the r„ for n > 2 are polynomi- 
als in q (with coefficients being rational in a, (5 and the 
four initial data). So far we have only found an induc- 
tive argument for this, but ideally Tn should be given by 
a determinant, as for linear q-special functions [9j. When 
the formula ( 11.51) for the solution to the autonomous recur- 
rence ( 11.41 is specialized to rational curves, then there are 
sequences where r„ satisfies a linear recurrence. However, 
for ( 13.21 ) it is possible to prove that this can only occur for 
the special values q = ±1, ±?. 

The recurrence relation ( 13.21) is invariant under gauge 
transformations r„ -^ AB"^ Tn-, while the transformation 
Tn -^ C" Tn mcaus that a can be rescaled to 1. Let us 
consider the particular sequence ioi a = (3 = 1 defined by 

g" Tn+lTn-l + (r^)^ 

Tn+2 = , T_2 = r_i = To = Ti = 1, 

Tn-2 

which is a non-autonomous version of the Somos 4 sequence 
( 11.21) . We would like to make the following observations: 



• 



• 



• 



For n > — 2, Tjj = Tn{q) is a polynomial in q of degree 
dn = nVl8 - n/6 + ^ sin(27rn/3). 

The coefficients of Tn{q) are all positive integers; the 
coefficient of g*^" is 2[^l 



The roots of Tn{q) appear to lie in a circle of maximal 
radius \/2 in the complex q plane, and coalesce to cer- 
tain special values, one of which is —1.3877957 to 7 
decimal places. 

The rigorous analysis of the properties of the qdPI tau- 



functions will be the subject of future wor 



^During the SIDE 6 meeting, Sasha Bobenko and Vassilios Papageorgiou informed me 
that this should be related to the Laurent phenomenon for certain recurrences recently 
studied by Fomin & Zelevinsky, arXiY. |math.RT/0311493| 
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